Replica symmetry breaking transition of the weakly anisotropic Heisenberg spin glass 

in magnetic fields 



m 
O 
O 

(N 
-i— > 
O 

O 

in 



S3 

c 

I 

C/3 



C3 



I 

c 

O 

o 



> 

00 

m 
o 

m 
O 

c3 



c 

o 
o 



X 
S3 



Daisuke Imagawa and Hikaru Kawamura 
Department of Earth and Space Science, Faculty of Science, Osaka University, Toyonaka 560-0043, Japan 

(Dated: February 2, 2008) 

The spin and the chirality orderings of the three-dimensional Heisenberg spin glass with the 
weak random anisotropy are studied under applied magnetic fields by equilibrium Monte Carlo 
simulations. A replica symmetry breaking transition occurs in the chiral sector accompanied by the 
simultaneous spin-glass order. The ordering behavior differs significantly from that of the Ising SG, 
despite the similarity in the global symmetry. Our observation is consistent with the spin-chirality 
decoupling-recoupling scenario of a spin-glass transition. 



Whether or not spin-glass (SG) magnet exhibits a ther- 
modynamic phase transition in applied magnetic fields 
has been a long-standing issue0. This issue is closely re- 
lated to the fundamental question of whether the SG or- 
dered state in zero field accompanies an ergodicity break- 
ing not directly related to the global symmetry of the 
Hamiltonian, i.e., the replica symmetry breaking (RSB). 
Most of numerical studies on SG have focused on the 
properties of the Ising SG0. Since the Ising SG pos- 
sesses no global symmetry in magnetic fields, the occur- 
rence of a phase transition in fields would necessarily 
mean the occurrence of RSB. Unfortunately, numerical 
simulations on the Ising SG have been unable to give a 
definitive answer concerning the existence of a SG tran- 
sition in magnetic fields 0, Q 

Experimentally, some evidence against an in-field tran- 
sition has been reported for the strongly anisotropic 
Ising-like SG Feo.sMno.sTiOsyJ. In fact, many of 
real SG materials are more or less Heisenberg-like 
rather than Ising-like, in the sense that the magnetic 
anisotropy is considerably weaker than the isotropic 
exchange interaction^]. Recent experiments on such 
weakly anisotropic Heisenberg-like SG suggested the oc- 
currence of an in-field SG transition[£j, in apparent con- 
trast to Ref.0. 

Meanwhile, via recent theoretical studies, it has be- 
come increasingly clear that the Heisenberg SG possesses 
an important physical ingredient which is totally absent 
in the Ising SG, i.e., the chiralityvL 0- In particu- 
lar, the chirality scenario of Ref.Jy, bj claims that the 
chirality is a hidden order parameter of the SG transi- 
tion of real Heisenberg-like SG magnets. According to 
this spin-chirality decoupling-recoupling scenario, in the 
fully isotropic Heisenberg SG, the spin and the chiral- 
ity, which are coupled at short length/ time scales, are 
eventually decoupled at long length/time scales, and the 
system exhibits a chiral-glass transition without the stan- 
dard SG order. In the weakly anisotropic Heisenberg 
SG, the Heisenberg spin is "recoupled" to the chirality 
at these long length/ time scales via the random mag- 
netic anisotropy. The SG order of the weakly anisotropic 
Heisenberg SG is then dictated by the ordering of the 



chirality. In zero field, some numerical support for such 
a scenario was reported [il Wi [ToL UlL fl^] . although some 
other groups claimed that the chiral-glass transition of 
the isotropic system already accompanied the SG order 

mm. 



In connection to experiments, in-field ordering prop- 
erties of the Heisenberg SG, both with and without the 
random magnetic anisotropy, are of obvious interest. The 
recent Monte Carlo (MC) simulation of the fully isotropic 
Heisenberg SG revealed that the chiral-glass transition 
corresponding to the breaking of the chiral Z2 symme- 
try occurred even under fields, while the observed tran- 
sition line had a striking resemblance to the so-called 
Gabay- Toulouse line observed experimentally^^]. Note, 
however, that the fully isotropic Heisenberg SG possesses 
Z2 x 5*0(2) symmetry even in magnetic fields, the chi- 
ral Z2 corresponding to a global spin reflection with 
respect to the plane containing the magnetic-field axis 
and the SO(2) to the global spin rotation around the 
magnetic- field axis. In the more realistic case of the 
weakly anisotropic Heisenberg SG, there no longer re- 
mains any global symmetry under fields. Hence, from 
symmetry, the situation is common with that of the Ising 
SG with the infinitely strong anisotropy. Meanwhile, the 
Heisenberg SG possesses the nontrivial chiral degrees of 
freedom, which are totally absent in the Ising SG. 

Under such circumstance, it would be highly interest- 
ing to examine the ordering properties of the weakly 
anisotropic three-dimensional (3D) Heisenberg SG in 
magnetic fields as compared with those of the Ising SG. 
Such an interest is further promoted by the apparently 
contradicting experimental observations on the Ising-like 
and the Heisenberg-like SGs0, Q. In the present Let- 
ter, we study by extensive MC simulations the spin-glass 
and the chiral-glass orderings of the weakly anisotropic 
Heisenberg SG in magnetic fields. 

The model considered is the anisotropic classical 
Heisenberg model on a 3D simple cubic lattice defined 
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by the Hamiltonian, 



A' 



(ij) 



(i) 

where Si = (Si X , Si y , Si Z ) is a three-component unit vec- 
tor, and H is the intensity of the magnetic field applied 
along the z direction. The isotropic nearest-neighbor ex- 
change coupling Jjj is assumed to take either the value J 
or — J with equal probability while the nearest-neighbor 
random exchange anisotropy 's (/i, v=x 1 y,z are 
spin-component indices) are assumed to be uniformly dis- 
tributed in the range \—D : D], where D is the typi- 
cal intensity of the anisotropy. We impose the relation 

T)^ v = = D u>1 

ij ji ij ' 

Simulations are performed for the anisotropy D/J = 
0.05 and for a variety of field intensities in the range 
Hj J = [0.02:5.0] , whereas the results shown below are for 
the field intensity H/J = 0.05 at which most extensive 
calculations are made. The lattice contains L 3 sites with 
L = 4, 6, 8, 10, 12, 16 and 20 with periodic boundary 
conditions. Sample average is taken over 40-600 indepen- 
dent bond realizations. To facilitate efficient thermaliza- 
tion, we combine the standard heat-bath method with 
the temperature-exchange technique . Following the 
procedure of Ref.^li care i s taken to be sure that the 
system is fully equilibrated. 

The local chirality, Xiy-i is defined at the i-th site and 
in the /i-th direction for three neighboring Heisenberg 
spins by the scalar, Xip = ^i+e^ ' x &i-e )' w here 
(/i = x,y,z) denotes a unit vector along the /i-th 
axis. By considering two independent systems ("repli- 
cas") described by the same Hamiltonian JJ), we define 
the overlaps of the chirality and of the spin by, 
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where (1) and (2) denote the replicas 1 and 2. In our 
simulations, the two replicas 1 and 2 are prepared by 
running two independent sequences of systems in par- 
allel with different spin initial conditions and different 
sequences of random numbers. 

We first compute the equilibrium time-correlation 
functions of the chirality and of the z-component of the 
Heisenberg spin (henceforth simply denoted as "spin" ), 
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where (• • ■) represents the thermal average and [• • •] the 
average over the bond disorder, while the "time" t is mea- 
sured in units of MC sweeps. Note that, in the above def- 
initions, the second terms [(g£)] and [(q^)] have been sub- 
tracted, which are nonzero even in the high-temperature 
phase in the L — > oo limit due to the absence of the 
global symmetry. This subtraction guarantees that both 



qx\t) and ql z> (t) decay to zero as t — > oo in the high- 
temperature phase. In the possible ordered phase, by 

(2) (2) 

contrast, both q^ (t) and q z (t) decay to zero if the 
ordered state does not accompany the RSB, but tend to 
finite positive values if the ordered state accompanies the 
RSB. The latter property arises because in the presence 
of RSB the t — > oo limits of the first terms are gener- 
ally greater than [(qt)] and [(<^)]. In computing the first 
terms, we perform simulations according to the standard 
heat-bath updating without the temperature-exchange 
procedure where the starting spin configuration at t = to 
is taken from the equilibrium spin configurations gen- 
erated in our temperature-exchange MC runs, while the 
seconds terms are evaluated in the temperature-exchange 
MC runs. 

The chiral and the spin time-correlation functions of 
the size L = 16 are shown on log-log plots in Figs. 
(a) and (b), respectively. A comparison with the L = 20 
data indicates that the data can be regarded as those 
of the bulk in the time range shown, since no appre- 
ciable size effect is discernible either in the first or the 
second terms of eqs. (4) and (5). As shown in the insets, 

(2) (2) 

both (t) and q z (t) exhibit either down-bending or 
up-bending behavior depending on whether the temper- 
ature is higher or lower than T ~ 0.21 (in units of J), 
while just at T ~ 0.21 straight-line behavior correspond- 
ing a power-law decay is observed. This clearly indicates 
that both the spin and the chirality exhibit a phase tran- 
sition at T = T g ~ 0.21 into the low-temperature phase 
where the replica symmetry is broken. A simultaneous 
ordering of the spin and the chirality seems quite natural 
in the presence of the random anisotropy. The estimated 
transition temperatur e T„ = 0.21 ± 0.02 is lower than T g 
in zero field, T g ~ 0.24 [lalT^. indicating that the applied 
field suppresses the chiral-glass (spin-glass) ordering for 
weak fields. This suppression of T g due to applied fields 
is consistent with the experimental observation 0. We 
also calculate the transverse (perpendicular to the ap- 
plied field) spin time-correlation function, observing the 

(2) (2) 

qualitatively similar behavior as q^ and q z (the data 
not shown here). 

We then perform a dynamical scaling analysis both for 

(2} f2) 

qx (t) and q z (t), with setting T g = 0.21. As shown 
in Fig. n](a), the chiral time-correlation function scales 
very well both above and below T g with /3qg = 0.9 and 
zcg^CG = 5.5, where flee, v cg and zcg refer to the 
[(ff?*)] (5) order parameter, the correlation- length and the dynam- 
ical chiral-glass exponents, respectively. The estimated 
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FIG. 1: Temporal decay of the equilibrium time-correlation 
functions of (a) the chirality defined by eq.(4), and (b) the 
spin defined by eq.(5), for D/J = 0.05 and H/J = 0.05. The 
lattice size is L — 16. The data at T = 0.14, those in the 
temperature range T = 0.18 ~ 0.24 with an interval of 0.01, 
and those at T = 0.26 and 0.28 are plotted (from above to 
below in the insets). The estimated transition temperature 
is T g ~ 0.21. Main panels represent the dynamical scaling 
plots, where the exponents are taken to be (a) (3cg = 0.9 and 
zcg^cg = 5.5, while (b) /3sg = 0.5 and zsGt'SG = 5.1. Slope 
of the straight dashed line is equal to 2/3/zv. Insets represent 
the raw data, where the data at T — T g are given in red. 



chiral-glass exponents are not far from the correspond- 
ing zero- field values [l^. By contrast, the spin time- 
correlation function scales not quite well for any choice of 
the fitting parameters: An apparent "best fit" is given in 
Fig. n (b) , whereas the quality of the fit is not improved 
for other choices. Such a poor scaling suggests that, 
for the spin, the asymptotic scaling regime has not been 
reached in the investigated time range and that the fitted 
values of (3sg and zsg^sg might not be true asymptotic 
values. The reason why the spin time-correlation exhibits 
a poorer scaling even in the regime where the chirality 



exhibits a nice scaling may naturally be understood from 
the spin-chirality decoupling-recoupling scenario: In this 
scenario, the chirality is always the order parameter of 
the transition, while the spin, decoupled from the chiral- 
ity in the absence of the anisotropy at long length/time 
scales (estimated to be r > 20 and t > 10 5 ), is recoupled 
to the chirality in the presence of the anisotropy at these 
long scales and eventually shows essentially the same or- 
dering behavior as the chirality, x^S. At shorter scales, 
by contrast, the spin and the chirality are trivially cou- 
pled by its definition, roughly being \ ~ S 3 , irrespective 
of the anisotropy. If so, an asymptotic critical behavior 
of the spin would be evident only at longer times t ^ 10 5 , 
which is beyond the time range of Fig. 2] 




FIG. 2: Temporal decay of the equilibrium time-correlation 
function of the spin, defined by eq.(5), of the ±J 3D Ising SG 
in a magnetic field of H/J — 0.05. The system size is L = 
16 averaged over 125 samples. The data in the temperature 
range T = 0.60 ~ 1.4 with an interval of 0.1 are plotted 
in the main panel from above to below. The data of the two 
lowest temperatures T = 0.60 and 0.65 are shown in the inset, 
together with the fitted straight lines. At any temperature 
studied, no up-bending behavior is observed. 

For comparison, we also calculate the spin time- 
correlation function JSJ for the 3D Ising SG with the ± J 
coupling for the field H/J = 0.05. The result for L = 16 
is shown in Fig. |2 Here, the behavior of qi 2 ^ (t) differs 
significantly from that of the weakly anisotropic Heisen- 
berg SG: Although the temperature range studied is as 
low as about 60% of the zero-field transition temperature 
T g (H = 0) ~ 1.1, which is expected to be deep in the or- 
dered state according to a tentative estimate of Refs. 
no clear up-bending behavior as observed in the weakly 
anisotropic Heisenberg SG is observed. Instead, q z (t) at 
lower temperatures persistently exhibits an almost linear 

(2) 

behavior. This is illustrated in the inset where the q\ (t) 
data at the two lowest temperatures are shown. In con- 
trast to the Heisenberg case, a comparison of the L = 16 
data with the L = 12 data indicates that some amount 
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of finite-size effect still remains in the second term of 

(2) 

eq.(5), and hence, in q\ '(t) itself. Nevertheless, the al- 
most linear behavior without any discernible up-bending 
tendency is robustly observed in common for L — 12 and 
L = 16, suggesting that this feature is a bulk property. 

In fact, a similar linear behavior was also observed in 
the 3D Ising SG in zero field^3> where the existence of a 
finite-temperature SG transition is established^. Unfor- 
tunately, we cannot tell from Fig. [^whether a true phase 
transition occurs or not in the 3D Ising SG in a field, since 
an apparently linear behavior could also arise when the 
correlation time simply exceeds the time window of the 
simulation, but stays finite. Yet, we may safely conclude 
here that the ordering behavior of the Ising SG differs 
significantly from that of the weakly anisotropic Heisen- 
berg SG, despite the similarity in their global symmetry 
properties. 

In order to get further insights into the nature of 
the RSB transition of the present model, we show in 
Fig. the chiral-overlap distribution function P(q x ) = 
[{S(q x — q' x ))] at T — 0.18, somewhat below T g . In addi- 
tion to the primary peak corresponding to q x = q x A > 0, 
which grows and sharpens with increasing L, there ap- 
pears the second peak at q x ~ 0, which also grows and 
sharpens with increasing L. The existence of two distinct 
peaks, both growing and sharpening with increasing L, 
is a clear indication of the occurrence of RSB. As ob- 
served in Ref.||, Pxilx) m zer0 exhibits a feature 
of one-step-like RSB, a central peak at q x = coexisting 
with the self-overlap peak at q x — q x A . The P x (q x ) ob- 
served here may be regarded as the in-field counterpart of 
the zero-field P x (q x ), with a feature of such one-step- like 
RSB. 

0.045 




FIG. 3: The chiral-overlap distribution function at a tem- 
perature T = 0.18 for D/J = 0.05 and H/J = 0.05. The 
transition temperature is T g ~ 0.21. 



In summary, by performing a large-scale equilib- 
rium MC simulation, we have shown that the weakly 



anisotropic 3D Heisenberg SG in magnetic fields exhibits 
a thermodynamic RSB transition in the chiral sector, 
which accompanies the simultaneous SG order. The 
ordering behavior of the weakly anisotropic Heisenberg 
SG differs significantly from that of the Ising SG de- 
spite the similarity in their global symmetry properties. 
Our observation is fully consistent with the spin-chirality 
decoupling-recoupling scenario, and might give a clue to 
resolve the apparent experimental discrepancy between 
the strongly anisotropic Ising-like SG and the weakly 
anisotropic Heisenberg-like SG. 

The numerical calculation was performed on the HI- 
TACHI SR8000 at the supercomputer system, ISSP, Uni- 
versity of Tokyo, and Pentium IV clustering machines in 
our laboratory. The authors are thankful to H. Yoshino 
and K. Hukushima for useful discussions. 
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